For Reference 


= 
-) 
io) 
pct 
DM 
‘com | 
se) 
SI 
= 
fe) 
oa 
ea 
Zz, 
ea 
ns 
< 
SI 
--) 
2) 
ke 
i) 
Z 


Ex IIBRIS 
UNINERSTCACIS 
AIPERCALNSIS 


= l= 
SEs || aaa 
| 
alt th a wy 
He ete 
Ww et* 
ug 


SDSS 

The University of Alberta | 
Printing Department 
Edmonton, Alberta 


ee wa 7 
4 A ie) ' 
AW Wy i] 
oO 4 ’ nk ut en 
a , 
i s v 


+ f 


me um mas 


" - ny wv + 


Digitized by the Internet Archive 
In 2024 with funding trom 
University of Alberta Library 


https://archive.org/details/Hertog1973 


ss  Gontieeana (alk gor 


_ a 
Ti Fe 


—_ ¥ oe 
a : eat SRC GRANTS mye 


_ Paumission Se haseky pene 
, AEGAN. Leena ae eure anes ey at 


_ aos to lend or Sel} poci-vinph eats al 
a = scientific re osm 


Dn Gah RN, Ceey ion Rm same ls One ALB BROT A 


RELEASE FORM 


NAME OF AUTHOR: William den Hertog 
eee Or Trine tos Excited State Wave Functions by 


Configuration Interaction 


YEAR THIS DEGREE GRANTED: 97S 


Permission is hereby granted to THE UNIVERSITY OF 


ALBERTA LIBRARY to reproduce single copies of this thesis 


and (to lend sor sell such copies forvprivate, scholarly, 


or scientific research purposes only. 

The author reserves other publication rights, anc 
neither the thesis nor extensive extracts from it may 
be printed or otherwise reproduced without the author' 


written permission. 


S 


Hes, Srp’ * 


Pee ee 


mL eSct Soa | : 


vad 
C7 
VAT iS4 4 = 


iivedasy TANS thee ill 


sity e@erded 0) andueteav’ 
vl A's ee hs | 
7 vy (Tie Ss sh Sel a Dew 


6 i] SEDAGT +} }J2meloa 


THE UNIVERSITY OF ALBERTA 


EXCITED STATE WAVE FUNCTIONS BY CONFIGURATION INTERACTION 


by 


WILLIAM DEN HERTOG 


A THESIS 


SUBMITTED TO THE FACULTY OF GRADUATE STUDIES AND RESEARCH 


IN PARTIAL FULFILMENT OF THE REQUIREMENTS FOR THE DEGREE 


OF MASTER “OF 7SCLENGE 


DEPARTMENT OF CHEMISTRY 


EDMONTON, ALBERTA 


FADO, sL973 


: ‘ = Vee lee 
ORAS oi, adore seaman 0 am min ‘9 
aaah axr Adt on Ba SNE 

a2471d2 20 egean 40 ae 


THE UNIVERSITY OF ALBERTA 


FACULTY OF GRADUATE STUDIES AND RESEARCH 


The undersigned certify that they have read, and 
recommend to the Faculty of Graduate Studies and 
Research, for acceptance, a thesis entitled 
EXCITED STATE WAVE FUNCTIONS BY CONFIGURATION INTERACTION 
submitted by William den Hertog in partial fulfilment of 


the requirements for the degree of Master of Science. 


WaKadad Ui vatobite SUAUORMD "x0 -THANOAA 


: 
: 1 
fas ,686s eved yodd tect? y2lesen bonpleteinw edt a 
: 
ten eaibuts- sysyhexP Io ytluess ody O38 Saommiooss an 
if ae fettizne wheats « ,etRntqeoos x03 ,fotesesh | 
a 


MOPTIASTUL WORTEMUNTANND ¥E BAOLWMIT SVAN ARATE GATT OK 


~ 


Within eheir irame: layup, 12 


At least those atoms whence 
TOMER OW obo L£are and send 


To shoot the sparks and scat 


Of the Nature of Things, 


ee rer en em nn an nd er Atm a me ce ae nF are 


x 


(trans letvon by William Bile 


Fee 
F- 
h- 


s by fire consumed ablaze, 


ie tr Raa 15 RIE PeAGeG) \ at ne) 


ABSTRACT 


The problem of obtaining wave functions and 
energy bounds is treated in this work within the context 
Ovecve cont laura. on 1nceraction (Cri method. Two dis-— 


Giuctemethods in obtaining the Cl function or, for the 


cr 


k state, ork with true energy Whee are investigated: 
(1) samul eous Minimization of the "energy error" 


~ 


2 A 
> and of ey / Ay 


Z eee Lier ; 
Wd | o> and (2) minimization of 


above by the Rayleigh-Ritz method. 


with the “local energy 
= <6, | (H - 
ae 
The first method is found not to give dependable 
numerical bounds to the deviation of one from iid anda 
new analysis of- the reason for this is discussed. Thus 
the major emphasis of the work is upon the use 


the Rayleigh-Ritz method in obtaining excited state wave 


in 


functions and upper bounds to energies. It is shown that 
the results thus obtained for two-electron ground states, 
ywery Closely satisfy the criteria of the first method. 
The group of states chosen for this model investi- 
gation is that of the lsns Ts states of helium. A systematic 
search for a useful CI basis has culminated in the 
Production Of ww Set Gf CI functions for schese states 
WLctrenerag ves within-0.003 to 0.005 2.u, of the corres-— 


ponding true state energies. Shortly after the completion 


of this work superficially similar results were published 
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but the basis used there introduces errors which are not 
: ul 
present in. the results reported here.. Thus .the lsns “S 
state functions and energies obtained in the present 
work are not only the most dependable yet obtained by the 
Rayleigh-Ritz method, but also show the promise of the 
method for obtaining excited state wave functions. 
The present work stimulated the development, by 
Dives mot lice ueativana | Soubical” analysis, whacheis 
similar to the well established natural orbital analysis 
but reflects the relative spatial concentration of the 
two electrons in the states treated in this work. The 
analysis is presented as used here for the first time. 
The expected decrease in importance of correlation and 
: . + 
approach of the inner electron function to a He 1s wave 
tunction upon consideration of progressively higher 
excited states is placed upon a firm numerical basis 
by the analysis. 
A principle is: established for constructing a CI 
; Bs Aone Te; his nae : 
basis toOrvocubiy. excited’, S state functions... “Such a basis 
must be formed from configuration functions spanning-a 
ae Les ee : ; tb ee a 
Maxiiial mark of the space or Lsns S state functions as 
well as from contiguration ftunctaons lying as much as 
s £ il ted ‘ r 
DOsSsivicwim ohe Sspece Of -2Zsns  S state Eunctions. The 
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Pe “Theory of Variational Metnodas 


ee Introduction 


It is well known that the time-independent, non- 


t 


relativistic Schroedinger equation of free atoms or 

ions with two or more electrons defies exact solution. 
im ehemwnecnod Of Conriguration interaceion, applied an 
Che work €O De 2eported, an approximations Lo the e@xactk 
SOlUEION coBpsaStarolt fa trial ftinctilon constructed on 

the principles of vector algebra and the calculus of 
variation. The function is set up as a linear expansion 
of a selection of basis functions expected to best span 
the space containing the exact solution aimed for. 

A Geunetronal of the triay “function wis chosem torserve 

as a criterion of goodness of the function. The latter 
is improved by optimizing the expansion coefficients 

and parameters of the basis functions to obtain a 
Stationary value for the functional as close as possiblie 


to the exact value known by experiment or theoretical 


In the subject of this thesis there are two func 


of main’ concern. With the Schroedinger equation written 


as 
HY, = WY 
k a 
ii wiaecky & = 0 refers to the ground Susate end k > <) peters 
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: th : ‘ : 
to the k excited state of a certain symmetry, and with 


a trial function to approximate ¥, denoted by 6 


k ke 


these functionals are defined by: 


k 
Ay. = <0, | (H W,)*|o,> 
They define: 
aes (E) -Wy) 
and Scere 


The definitions imply that Ey and Ay > 0 as OD eas Yee 


The thesis reports on the study of two different 
Vartational methods. to construct excited State wave functions. 


The methods are based on: 


Poe ene Minin ZatLon Or ey 


ae 
é : A Se ie fa 
2. The simultaneous minimization of eA, and Ave 


I.2 The Minimization of «©, 


The configuration interaction method, applied in this 
work, is the Rayleigh-Ritz variational procedure in which 
thedbasiseftunctions are «symmetry sadapteditwo-electron 


configurations. The method, mainly appilved to. qround 
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states, owes its usefulness to the fact that the addition 
of configurations in the trial function is equivalent to 
approaching the use of a complete basis set. 

For the application to excited states the MacDonald 
theorem is of central importance (1). Applied to the 

: : : ; ey ; 
Ci method it establishes that the m lowest eigenvalue 
obtained from any given basis set of N symmetry adapted 
Oh aren ae rea th 

configurations 1S an upper bound of the exact m lowest 
energy of the states of the symmetry under consideration. 

Clearly the theorem immediately provides a variational 
criterion, namely the minimization of ce, = EL, - W, or Of 

ass 
Dm _ rT - te an ; 4 ie 1 1 x 
Bh 1£ the exact energy of the k excited state is not known. 
‘ 

Yet, while the theorem was established as long as some 
thirty years ago, and since has béen occasionally referred 
COwinecie texts “and di teraturer (2.3,;4,5)> ne application 
to excited states was found in the literature at the time 
that the work for the thesis was undertaken. 
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An explanation of this unexplored opportunity is 
Giplricule to Gave. .ltrmay lie in. the tacusthat-in the 
large amount of ground state work done, the linear 
variations to obtain the lowest Ea value at the same time 
would yield a set of BE. values invariably too high to be 
of any use as upper bounds. For example, the basis set 


ised by Weiss (4) generates for helium in the same 


Minimization a ground state energy of -2.90270 a.u. (accurate 
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Sub Oe40ed. te tO, the eGxact energy Of  -2elt4o0> alu. Of the 

ban Wp atoge aS excited state. This may have led to a general 
impression that the MacDonald theorem was not of much practical 
value, although in that case it must have been overlooked 
thatmatbasis «set.of configurations good. for: a ground 

state must, of course, be expected to be a poor one for 

excited states. 

Thevaudbetous=choice,of.a particular basis ~setsis 
crucial if the MacDenaid theorem is to be applied to 
excited states. Indeed this consideration has been the 
starting point of the work to be reported below. Several 
aspects have been explored such as different types of 
orbital functions in the basis configurations, the size 
of the basis sets in relation to the problem of 
convergence towards the exact energy, the role of the 
linear (expansion coefficients) and non-linear (orbital 
exponents) variational parameters, and the selection of 
the configurations. 
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These studies have led to the computation of a set 
1 ny : plas oh 7 
of ~S states of helium with energies sufficiently close 


to the exact ones to demonstrate the strong potential of 


the Rayieigh-Ritz) variational method,for excited states. 
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I.3 The Simultaneous Minimization of ef/h and A 


bi ea BEA e 3 


A method which uses the functionals Ay and Ee 
as variational criteria has recently been proposed 
by Messmer (6) and modified by Choi, Lebeda and Messmer (7). 
The method is presented as a better alternative to 
minimizing either Ey, OF A, alone. It had been shown 
(i332) that in the latter minimizations wave functions 
are obtained which have minimum errors in certain regions 
of configuration space, rather than over the whole of 


configuration space. On the one hand in minimizing Ay 
the errors in the wave function associated with the 
closeness of the electrons will be minimized; in this 
connection n is called the "local energy error" (13), 
-representing with reference to the closeness of the 
electrons, “Short range errors.” On the other hand, 
"long range errors" are minimized by minimizing Ens 
the “energy error" (13,14). In this case the wave functions 
obtained should be suitable for describing properties 
such as the radial density or diamagnetic susceptibility 
baal long range errors are important. 

The basic premise of Messmer's proposal is that 
the best approximation to the exact wave function hag is the 
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expression above illustrates how simultaneous minimization 
of afi and Ay should diminish the deviation coefficient 
as well as narrow its range, thereby yielding a trial 
fune tion Ds with maximum ay. It is evident that the 
usefulness of the method then depends entirely on the 
goodness of the proposed bounds. Even if a good oF would 
be obtained, its goodness can only be established if the 
bounds enclose a narrow range for ‘the corresponding 25 
ioe it Sepoints which requires closer consideration. 
Some work had been done which strictly as numerical 
evidence has thrown serious doubt on the goodness of the 
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example of the method is given in an application to the 
harmonic oscillator. Upper and lower bounds are 
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Ti. # Computational Precedures 


The Hamiltonian used here is the non-relativistic, 
spin-independent, infinite-nuclear-mass Schroedinger 


Hamiltonian for two electrons: 
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Separate computer programs were developed for the 
Minimization work with the Laguerre and Slater type orbitals. 
The correctness of both programs was confirmed by com~ 
plete duplication of published Helium ground states on 
bases of twenty configurations, respectively by Holoien 
(9) and Weiss (10). Ail computations were carried out 
in double precision on the IBM 360/67 computer of the 
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Choi et al (7) originally proposed the simultaneous 
Minimization of ec / hy and Ny, in the context of self- 
consistent-field theory for excited states because then 
ey alone is not an adequate criterion for obtaining a 
good wave function... In “the present. work with configuration 
interaction, the MacDonald theorem does permit the use 
Cr e, as criterion since there exist well-defined lower 
bounds to each approximate state energy. It then becomes 
of interest to enquire whether there are eee discrep- 


ancies in the nature of the minimization process usin 
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the ey /A k and Ay versus Ey minimization methods. 
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COntne. Ground states ot He ,.dil -andeaBe ©. "Or comparison 
these were all carried out on the same configuration 
bacd seme Tha si) basis -aias set, often, coniagurations,, is basis 
VI listed in Table 10 and is expressed in the associated 
Laguerre functions described in Chapter II. 


* further details of the basis are given in Chapter IV. 
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The results are listed in Tables 5 and 6 and plotted 
in Frguress 374, "and 52 In rigures @3 (a)-64(a)s and 5 (a) 
each piotted point on the curves labelled EG represents 
the minimized 6 obtained by linear variation at the 
indi caAccadmvalue of na To each EG value belongs a trial 
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Thus the EG minimization produces three curves for a 
ground state: the Ee, curve obtained by variation and the 
- ee Fad 
A-and e« /A-curves calculated from the trial functions 
corresponding to the ¢« _-curve. Similarly, “as Shown in 
Figures 3s (ie, 74b)s,; and S(t), the AG Minimization leads 
to three curves: the A curve by variation, the ex 
Ne es ‘ bin 
Curves and. « /a—-curve by calculation Erom’ the triad 
functions Gerresponding to the A-curve. -For comparison 
all minima are given in Table 7. 
The following observations are made: 
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For He the lowest p and Ay Soa nende, at iy e= 2.325% 
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Ao and he werescaloilated. — ~ltups. secnminetpesrabplae 


thaeopne «-lowerethe value .of a4 the lower are the values 


ore : : 
of Ao and Ef be As in. the preceding case the criteria 
of lowest Ey and lowest cf / by run parallel. By either 
erilerron the four trial functions wouldsbe arnuangedian 


the same order of goodness. 


Eb. CONC Lustons 


The above analysis leads to the following decisions. 


First, the minimization of A, alone is rejected in favour 


k 


fe EL Mie nNinuzZactOn .w. Lhe reasons aveathat the isimin tancous 


decrease of EL and Ay is best approached in the Ey 


Minimization lobservations, (i) and..(11). above), and that 


the lowest values for oe and ct Fees are obtained by e 


Minimization, observation (ive above). 


k 


Secondly, in the EL minimization the calculations 


of Ay and ie can be omitted since they do not 


significantly add to the information provided by the ey 
Minimization [observation(iii) above]. This has also 
the great advantage of cost reduction since the 
Galcibation ,of dy requires far more computer time than 
Goat. ot Exe 

In conclusion: The Rayleigh-Ritz variational method 


is found to very closely satisfy Messmer's criterion 


that en hins and dy should decrease Simultaneously. This 
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being so there appears to be no gain in using the 
methods involving Ay as an alternative to minimization 
of Ex in context where the MacDonald theorem makes a 


minimization a viable technique. 
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III.4 The Upper Bound 4, /(W} - Wy)” 


+ + 
the work, ongthe ground states Gf He, Li’, and Be “ 


reported above produced some data which are pertinent 
to the upper bound A, / (wt - W,)*. 

It was pointed out in Chapter I that the weakness 
of the bound follows necessarily from the spacings in 
the atomic energy spectrum. This was illustrated with a 
hypothetical energy spectrum. Some supporting numerical 
results can now be presented. These are obtained from 
the eo minimizations, reported above, on the same 
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ete 3 
basis for He, Li , and oy (fable 7). FErom the trial 


functions belonging to the Eg 8 the values of AG are 


calculated, and since the accurate values of |W, - Wo 


ak 
are known (17), the upper bounds Ao Ww, - wo)? can be 


calculated. Mhe results are Twisted in Table 39. 


ore - 9 
me defined by A, = <o, | (H Wy) bey 


Ay + 0 as oy > Yee Bearing this in mind one would conclude 


from the A. walues in*Table 9 that, for example, ec of 


Since A 


He with A, =e. O2 72 shemld bewa mMOregaccurate trial 
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upper bounds would lead to the opposite conclusion. 


These pounds, petngs:0:..4 283 ston He jancga0. 0708 for Bee 


imply that the squared deviation coefficient as k Ox 
v 


oy for, He could bes six times thatsor a for Boe 


The reason for this reversal lies in the values of 
++ 
|W. - We! as can be seen in the table. Since Be 


has the largest |W = We! value and since the upper bound 


dt 
; : 2 ++ 
is inversely proportional to (We = We the oy of Be 


Lig. LOsLes fe Val vemand 
O 


Cu 


would appear the poorest accor 


, 2 
WieeS / — WW x 2). 
the best by its Af WW, ee value 
These results therefore clearly illustrate that 
AL / (Wy - W im cannot be expected to be a dependable 


upper bound. 


IV. el ne vapplicattion ofere 


Rayleigh-Ritz Variational Method 


LV oie introduction 


As discussed in the preceding chapter a comparative 
study of Ey and Ay minimizations decided in favor of 
applying the Rayleigh-Ritz variational method for the 
calculation of excited states. In this chapter the search 
for excited state wave functions based on this method 
will be reported. In the course of this study a number 
of bases were generated which, for convenient reference, 
are collected in Table 10. These bases represent a 
gradual improvement in their ability te generate trial 
piunctions with venergies close tomthe true energies. 

The study culminates in finding a basis which firmly 
establishes the great promise of the Rayleigh-Ritz method 
for excited state wave functions, 

Sportly after the conclusion of this work a 
publication by Wang and Wed une td appeared in the liter- 
ature «(15) dealing with excited state wave functions of 
helium obtained by the Rayleigh-Ritz variational method. 
This paper 18, thee tirst to deal withethevappii cation af 
thes “apprcacth in the literature, and” therefore merits a 
Separate; discussion below. “In this ‘connection one 
particular basis set investigated in this work is of 
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special importance since it shows flaws which help to 
explain inaccuracies in Wang and Weinhold's excited 


state energies. 


Ve The development of basis sets. 


the -Liret concern in this study was for obtain some 
idea Of the results one could get with the Rayleigh- 
Ritz variational method. For this purpose the choice 
of the bases could be arbitrary and fell on a basis which 
had been published by Holoien (9) as being a good one for 
the Ground Stace OPehelinn.—Tinsebacis NO. I toain 
Table 10) is expressed in the Laguerre functions described 
in Chapter II. With the orbital exponent optimized to 
2.222 the basis generates a ground state energy of 
-2.90123 a.u. which compares with Pekeris' (12) highly 
accUrace Valueror mas. 9037/2 a.u% 

Although this basis cannot be expected to give 
good approximations to the-excited states in the Ey 
minimization to be explored, the results are still of 
interest for two reasons. 

First, they illustrate the, type of minimization 
curves generally obtained in this work. As shown in 
Figure 6 a setiost Ey values is obtained by linear 


Variation for a specific value of «the orbital exponent 
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bounds Ey values were obtained. The minimum for each 
Ey was then. established by non-linear variation, i.e. , 
the Variavion of 7. As shown in the figure, a mini- 
Mizationecurvesis thus obtained for eacheE. value. 
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The minamarot the curves are listed in Table: Ul as En 


pairs. As to be expected the minimized EL values are 
useless as,upper bounds.,. All of, them, except Ey and 
lie evén above the series limit of -2.0 a.u. of the 
Singly excited I. states of helium. 

Secondly, the results illustrate a viewpoint stated 
earlier by the author. This was that generally in work 
on ground states the Rayleigh-Ritz method would yield 
upper bounds too high to encourage attempts to improve 
these for excited states. A case in point is provided 
by the upper bounds obtained with the minimization of 
Es ace =—s42.222. They are listed) ingthe fourth ‘colsmnp 
of Table 11. These Ey values are indeed too high to 
even remotely suggest a possible use as upper bounds. 
On the othereshand, it,is,worth nokingaschat thegnen-— 


linear variation does bring about an impressive lowering 


e) 


£ the upper bounds just mentioned (Table 11, columns 
Peano.) s eloaseats Nowever, bY ateeligetil lan taom 
sufficient, as evident from a comparison with the exact 
Values W. given in the last column of the table. 


Further improvement is therefore to be sought 
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in changing the basis itself, by varying and optimizing 
its composition and length. To study the composition 
with fixed length two sets each containing 10 con- 
figurations were selected from Holoien's basis just 
ouscuscseda(No..llt in Table 10)., The onesselection is 
bases No.1 which is simply the first 10) radial terme 
Oreoaals NO. Lit. Thesother selection 1s Dasis No. Ad. 
This set was selected by analysis of the trial function 
for the ground state obtained from basis No. III at 
lesececec. aerem this function *wererschosenm the ten 
configurations with the highest expansion coefficients, 
that is, the ten energetically dominant configurations. 
The results obtained from these sets aoe given 


£f interest specifi- 


e) 


gnevable l2tand;Figure 7. .They-are 
‘cally with respect to the most serious problem in CI 
calculations, namely the selection of those configurations 
which lead to a most rapid convergence toward a specific 
state. There is no general rule available for choosing 
Such configurations (3). Figure 7 shows the supériority 
Of Dasis No. VI over basis No. I in tezmms of con- 

vergence towards the true ground state energy. The se- 
lechion ineset No. Vi of the angular terms (2p), 
erates and isa)" in favor of radial terms leads here to 


a substantial improvement of the convergence. 
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Table 12. Ground state energies and eigen vectors; 


Basis. Setoel andevie 


Basis Set. y Basis Seta 
B92 8c 705974 Eu sre 899700) 
Oo O 
fee Deeg (0 n = 2.325 
Cio caueue? ere) ers 
2 2 
0.9612 mee 0.9029 se 
B25 46 lees 0.3940 Ve os 
0.1005 oe sty, eee 
| 2 
20.0322 be 0.0631 (2p) 
~0.0080 Pere 0.0304 vone 
: 
0.0052 ene 0.0272 peor 
~0.0049 heey 0.0199 eee 
0.0034 bale ONO 37 lees 
2 2 


#020026 (4s) 0.0127 (3p) 


(a Odo 3s4s Of ow 2s4s 


n* gh sap 
“reese = 


2 


“Ca s) 

ital 
‘aol 
ro 


oF a ae 


; 7 


» 


i ai ee 


a aa 


a 


; : a “ Ls ie - ha au 


{Lagry ° 
ae ee es 
o-8in0- 
010.9 
Ce 40- 
CRlo.i 


oon 


-2.8785 800 


from basis 
No 1 


-2.8786 006 sa al lt aden thease baat oe Sele 


cel 


-2.8997 500 : 
from basis 
No VI 


“O79 # £00 Fo 


\ 
* 
wy, 
| 


Mec ace, eee ant med oe ee ee ee 
1.94 7.06 


q 


PIgGure: j. bo iin ZaAtaon, fer He. 


5A’ 


important, it is the sequencing of these which is 

crucial to optimize convergence. In the absence of 
physical considerations to choose a sequence, the results 
With sets No. I and VE suggest an alternative. If a 
truncated basis is to be developed by trial and error 
then after each trial the configurations dominant in 

the trial function should be retained and the 

remaining configurations should be replaced by untried 
Ones corm the next trial. 

Ali trials, including those justedescribed, with 
orthonormal bases, that is, with sets expressed in 
Laguerre orbitals and one orbital exponents, are summarized 
in Table 13. Some of the results are noteworthy. 

The importance of angular correlation in the 
ground state was illustrated above with the two ten- 
term bases No. I and VI. Further confirmation is given 
by inemtact, Shown in Vable 13, that basis No. Vi with 
10 terms gives a ground state energy practically as good 
as basis No. IV with 20 terms. Interestingly, the eo 


term in basis No. VI is as effective as the eleven 
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(4s)~, 1is5s, 285s, 385s, 485s, (5s) , 2p3p, 2p4p, and 

3p4p. 


On the. other hand, for the excited states, the 


radial terms in the bases are more important. This is 
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seen in Table 13 by comparing Er Eve and E, obtained 
ErOmebases NOy DyPeand=iv. “The change =rom basiseNow rit 
to basis No. IV is equivalent to the replacement of the 
five angular configurations (4p) 7, (cane 3d4d, aan 
and (AEy DYwrne*radiakw tems lsusy "2s5se = ss5S, 4556" 
and (Sst. The result is that EL improves from -2.1336 
COere 420 *aeueeand EY Pron-—1 #9286 stor. 0245 "ae 

The improvement is, however, at the expense of the 
ground state: ES Changes ™ rom 259012 Sto 82750996 eae, 
consistent with the importance of angular correlation 
in the ground state discussed above. 

To follow the role of radial terms further a 
basis was then chosen in which angular correlation was 
left out entirely. This is basis No. V: Cis) aan 
Wsetsy iiscs, «s/f Uslss 2" Apartstrom-the Limst een nie 
consists of only singly excited state configurations of 
thestorm isis. This choice is seen in) Tabie-13 to 


“bring a drastic improvement which is the greater the 


higher the excited state: Bae Eo, fae Anas scotained 


NO 


with basis No. tll improve respectively by 0.01, 0.02, 
Cio paNe ic oO, Ghenkl is 

TWhesbases, aaiscussed Up to this point alls contain 
the restriction that the same variational value was 
assigned to all orbital exponents of one basis. Although 


further improvement would result from lengthening this 
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type of basis, greater efficiency should be obtainable 
by lifting the restriction to a common orbital exponent. 
Thus different values were chosen for the orbital 
exponents of basis No. III. 

For the Slee configuration the orbital exponents 
nee Zao ala Ny = 1.19 were chosen. These are the 
Optimized EGxpONnents in Hekert S varvation function (3,11) 
recognizing the radial correlation between the electrons 
Tim vole oCLoOuUnGs otare Ol Weimun. Fon =the: (isis) =con= 
figurations (n # 1) all 1s orbitals were assigned on 220 


and all ns orbitals n 1/n. This basis is con- 
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veniently called the 1/n basis. With a configuration 


r ; "ae 
represented by (n gVL mys) the basis can be written as: 
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The: 1 /mybasiomrs tabsos listed -injTable, > efor, later 


All preceding bases were orthogonal since the 
associatea Laguerre functions used as basis functions 
are orthogonal if they have the same orbital exponent 
(Chapter il). Since this advantagesot orlnogonality was 
Toot tnethe I/n set the Haguerre functions were meplaced 
by the simples Sliacerstype orbitals (Chapter 911): 


Since in the l/n basis all oxrbutal exponents are 
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faxedsonly a linear variation can be applied. The 
resulting Ey Values are listed in=Table 314 and there 
compared with accurate values Wye lt is clear thatethe 
1f/n basis has led to a dramatic improvement. 

The upper bounds agree with the Pekeris values 
LOuwitiinethestnird rounded off decima errom the; 4th 
me state upward. To the best of the author's know- 
ledge this is the first time that such close upper bounds 
were obtained by the CI method. 

The conclusion is that these results firmly 
establish the great potential of the Rayleigh-Ritz 
variational method for excited states of two electron 
systems. The results suggest also that it will be 
worthwhile to investigate the application of the method 


to excited states of systems with three and more electrons. 


LVe2° thee basis Versuce elem sl dmDasts 


Under this heading a digression must now be made 
betore the results obtained with the 1/n basis can be 
discussed in view of Wang and Weinhold's publication 
mentioned in the introduction. It was pointed cut that 
the b/nedasis permits only linear variation, Since all 
EPerorbiteal exponents have fixed valuess A possible 
modification which permits some non-linear variation 


is the use of a common variable orbital exponent instead 
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Thus, the 1/n basis was modified by replacing only 
the orbital exponent Ny = l/n by one variational value 
No: This new basis is called the No basvcw his 
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l/n basis: 
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The variation of No of the new basis yielded minimization 


curves for each Ey in the same way as illustrated earlier 
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States for k > 1. The upper bounds obtained with the 


Oe0Gu basis) are listed 1n Table 14 with those of the i/n 


One notes immediately that the upper bounds from 
the 0.08 basis are even lower than those from the 1/n 
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values of We. This Ey value is obviously not an upper 
bound and seems thereby to contradict the MacDonald 
theorem. But as mentioned in Chapter I this theorem is 
based son a wigorous prootethat: al. E.'S obtained by the 
Ravel gn Ritz variation must be upper bounds: to the 
corresponding Wi S- 

An explanation for this apparent contradiction 
was therefore sought in examining the accuracy of the 
computations involved. In the computational procedures, 
as described in Chapter Il, three general stages can 
be distinguished. 

i) the computations of the basic integrals in the 

Setting up of the Cl matrix. 

ti) the canonical transformation to orthogonalize the 
matrix eigen value problem [Chapter II, equations 

GB era ticla kes) 6 
naar} the diagonalization of the transformed CI matrix. 

Of these three the first and last could be ruled 
out as sources of large errors. In stage (i), in the 
calculation of the integrals no differencing errors 
occur since the expansions involved contain only 
addantions. “in stage \(i11) the classic Jacobi dragonala— 
gation is applied and is well known to be very stable (22). 


Errors could, however, be expected in the canonical 
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replacing theevalvesel/n and 1/n' assigned to 4 in’ the 
eEbitel functions (ns) and n's) by one constant such as 
0.08 will cause an increase of overlap between (ns) 
anda(nuis)ey Borckexample, in the orbital notation of 


Chapter iii: 
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Such an increase will affect the elements of the metric 
S of the configuration basis, and may cause errors in 
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the canonical transformation of S according to: 
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An orthogonalization check was therefore applied by 
scanning the off-diagonal and diagonal elements of I 

for deviations from respectively 0 and 1. The magnitude 
of the deviations should directly reflect the accuracy 
of the canonical orthogonalization. The following 


Maximum deviations in I from the true identity matrix 


were found: 
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DCs clear from these results that thescnorce 
One. Uo LOT nN» in the No basis causes an impermissible 
traccuracy ii tie OLenogonal 1 7zations sinestallure of Eg 
as an upper bound must therefore be ascribed to 
jazmuretrons?: of the computational precision. The 0.08 
basis must therefore be rejected. 

Another orthogonalization check, of importance 
for the discussion in the next section of Wang and Wein- 
hold's work, is obtainable by investigating the ortho- 
gonalization of the orbital basis as distinct from the 
configuration basis. It was found that the Schmidt 
orthogonalized orbital functions yielded an overlap 


matrix which showed the following maximum deviations 


from a true identity matrix: 
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Again the 0.08 basis Clearly Causes (Serious 
computational errors. In anticipation of the discussion 
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is noted here. An orthogonalized orbital function of 
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the 0.08 basis is an expansion in functions, most of which 
haven che esanetorp tra ltextonentis 1.6e., n = 0.08 -, 
whereas in the 1/n ‘basis’ an orthogonalized-orbital 


function is an expansion in functions with n defined 


differently for each function. 


TVA Comparison with Wang and Weinhold's Excited States 
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Wang and Weinhold's work introduced at the be- 
ginning of this chapter can now be compared with the 
results with the 1/n basis. For direct comparison their 
results are listed with the present results in Table 14. 
Wang and Weinhold give the numbers in parentheses as 
EhewEmiris (Of terror in the fina twee meciumals:. Seis 

of interest to note that their average Ey values are call 
below the Pekeris values, exactly like the Ea obtained 
from ithe 0.08 basis. The following discusston will 

show that these errors all have the same cause. 
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Their basis in terms of configurations and orbital ex- 


ponents can then be*written as: 
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and is quite similar to the 1/n basis which for comparison 


is repeated: 
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Wang and Weinhold comment on the rather large 
errors given in the parentheses in the table, which 
they obtain white carrying out the’ calculations in double 
precision on an IBM 360/67, exactly as in the present 
work.” They report that monitoring of the accumulation 
cf errors shows up rather large differencing errors in 
the repulsion integrals between the Laguerre functions. 
As they state these errors can occur wnen, for large 
quantum numbers n, the oscillatory behavior of the wave 
Concer ots becomes more and more pronounced. 

With nessect to the cause of these errore “a 
similarity between Wang and Weinhold's basis and our 
0208 basis willenow be discussed. The hydrogenic eigen 
; /n 
functions can be expanded in STO's (16). Let [n,2,m}\% ~ 
represent the hydrogenic function with quantum numbers 
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represent a STO in a similar way. The expansion has 


thesfornm: 
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NOtetnatathon orbitals exponents of athe STO! sia. se 1 /neand 
not 1/k. Thus the expansion generates (n-2) STO's with 
the same orbital exponent 1/n. .For example the orbital 
25S 2S een fact.a. summation .of 25. STO tic eau le waist hee ele 5. 
It is this abundance of orbitals with the same n which 
caused the errors in the 0.08 basis, and is evidently 
also causing the errors in Wang and Weinhold's results. 

It can therefore be concluded that the-choice of 
Slater type orbitals for the 1/n basis represents a 
substantial improvement over Wang and Weinhold's 


enolee of hydrogenic functions for therr basas. 
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Ve seine Cresunctions tor the 1, States of Helium 


Vi ime Introduction 


In the preceding chapter the Ey values were 
reported obtained with the Rayleigh-Ritz method applied 
to the 1/n basis. It remains now to report the corres- 
ponding Ci=functions for the poe states of He. 


There are different ways to present a CI function 


a The most direct one is as: 
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since the C31. '8 are obtained as eigen vectors in the 
diagonalization) of the CI matrix. 

An alternative is the well known representation 
igietatural “orbvtals {(NO'S) “introduced by Lowden (18). 
TheenactivaLl"orbitals x4 gare linear expansions in the one- 

< 
@reerrvon basis*functions, STO*s in the present work. 
These expansions are obtained by a transformation such 
that the superposition of configurations constructed 
from the natural orbitals has the most rapid convergence. 
Further, in two-electron systems, the configurations are 
all’ doubly o¢cupied orbitals. Thus the NO form of the CI 


Furection “is 
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riginal STO basis to the NO basis has been derived 


with a 


by Lowden (18) in the context of the theory of density 
matrices. An alternative derivation (19) is given below 
in section Ve2, This is of importance, with respect to 
Scce1onav. agin waich Yrational orbitatstswill begintro= 
duced and derived. 

The NO presentation for two-electron systems has 
been widely adopted because it ensures the most rapid 
convergence and also because the configurations being 
doubly occupied orbitals provide a description most 
suitable for ground states. But for singly excited states 
an expansion exclusively in doubly occupied orbitals does 
not seem. to be the most "natural" or rather "rational" 
expression. Indeed from a chemical rather than a mathematical 
point of view one would prefer a rapidly converging 
expansion of a such that cne configuration would be 
dominant and have singly occupied distinct orbitals 


R R 
X1k and Xo% for the inner and outer electron. The 
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expansion would then have the form: 
R R R * 
a as . ay 
% (1,2) rece’ aie [x (4) X (2) ] + e e ° 3) 


On the basis of these considerations Birss (19) 


has introduced the concept of "rational orbitals" 
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The square brackets will from here on refer to the symmetry 


adjusted form of the configuration. 
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(RO's). The development of this concept is given below, 
in section V.3, followed by an RO presentation of the 


CI functions obtained from the 1/n basis. 
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The natural orbital expansion can be obtained by 
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with a the "fundamental matrix" of Lowdin (1B) 


To obtain natural orbitals S is maximized subject 


to the normalization constraint 
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Using the method of undetermined multipliers (25) with 


Multiplier <«, One has 
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five equation 1s adentical to that obtained by Lowden (18) 
showing that the maximum overlap criterion does obtain 


the natural orbitals for the two-electron problem. 
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V..3 "Rational" Orbitals for Helium 


Consider the case where the natural orbital 


Representation Oied Cl functionris 
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This appears to imply doubly occupied orbitals. 


Indeed if lela a reasonable orbital approximation 
to the state is just Ix X71 So that, the’ naturale orbitals 
are rational in this sense. If, however, at and aly 

are of the same order of magnitude, either single 


determinant is a very poor approximation to the state. 


Just such a situation arises when the best orbital 
approximation to the state is 
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This can be seen by doing a natural orbital analysis 


on this function. The fundamental matrix is {assuming 
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Since (without loss in generality) the XG and A can 
be assumed to form orthonormal sets in themselves, the 


matrix R is orthogonal and can be written as, 
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Where does the natural orbital analysis “go wrong", 
iieenierscnse that Lt gives a function which Ts *practically 
Correct but analytically obscure? 

By the very act of diagonalizing the fundamental 
matrix one forces the expression of a CI function into 
a series of doubly occupied orbital configurations. 
What was done in the previous section was to "dis- 
dtagonalize" the matrix, using as criterion: the requires 


ment to maximize the off-diagonal element. To accomplish 
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the goal of obtaining a function which is analytically 
clear while containing some minimal number of terms, 
one wishes to transform the fundamental matrix so that 
it is zero everywhere except for a minimal block 
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off-diagonal element in that block as large as possible 
This procedure is ill-defined because it does 
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natural orbital analysis does specify a way to a 
minimal block (albeit diagonal): diagonalize the 
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Select ayo as the element to be maximized by 
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are presented in Table 17. The advantage of the RO 
presentation over the NO presentation is evident from 
the table: in all excited state functions the first 
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V.5 The Rational Orbital Functions 


In this section the expansion coefficients of 
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the highest coefficient would change from b._ to b. 
im in 
with n > m, this change would correspond to an increase 
of the average distance of the electron from the nucleus. 
In the light of above analysis the RO's can now 
be examined. They are presented in Table 18 and 19 
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Table 20 Characterization of Rational Orbitals ae 


k 
Dominant Orbitals Fe 
a pba inner electron outer electron 
au 
(lsns)'S COGLE. Ler Fe Crear ie sk Os kig F. 
1 0.983/0577 1 0.984 ih 
2 0.99906525 ‘i Ov SS 4 
3 On 9907 71 1 1. 1 0.924 5 
4 0.99991128 ue 0.949 6 
3S 0.99995747 1 0.927 i 
6 0.99997677 5 -0.877 8 
7 0.99998610 a On 12 9 
8 0.99999109 1 Oroo 10 
9 0.99999396 ne 20.65 L2 
10 0.99999571 il -0.741 3 
ime 0.99999682 al: 0.810 14 
iT? 0.99999753 iL 0.866 15 
3 0.99999793 2 0.910 16 
14 0.99999788 uh eons. i ey 
nes; 0.99894875 iL ? ? 


tiat the RO presentation’ of §theacz £finetion clearly 
reflects the relative spatial orientation of the 
electrons. 

As for the outer electron, there is, in Table 20% 


clearly a progression towards higher F, with higher 


k 
SyeleecQustate. By the analysis given above) of ‘the 
Schmidt-orthogonalization this means that RO describes 

the outer electron as being on the average further away 
from the nucleus as nin 1sns Is increases. 

In conclusion, the RO presentation of a two- 
electron CI function is in accordance with the chemist's 
picture of the inner electron being in a 1s He -like 
orbital. Also, as the state energy increases, the RO 
fopethe outer ¢@lectron indicates that the electron is 
progressively further away on the average from the nucleus, 


and reflects the corresponding decrease in correlation 


energy. 
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2 
Vive  ethe a 2s) 1, State of Helium 


VI.1. Introduction 


Compton and Boyce (27) observed an emission line 
in the vacuum ultraviolet, a corona line, with wave 


2 =] a3 aes 
number 279,715 cm ~~. Kruger identified this line with 
“hipaa ape rie A el Sa : : 
Chew transition (25) S = ls2s: Son the basis that, 
firstly, the line never occurred in helium spectra 


when another line, identified with the transition 


Reale 1 


CE is2s ~S-was absent and that, secondly, the 


f 


{ 


Hine fell within .a calculated 20,000 cm interval. 
No details of the calculations were given. Bransden 
and Dalgains (29) calculated on the basis of pertur- 


bation theory that the line may correspond to one of 
3 5 3 


THeecrtanistcions 2520, .P - ls4s SS, 2s2p > -P-— is5s Ss 
co 
+ BS) 3 
encecseo P= Is4q dD. 
24 * il ag ers 
Lieiwie. (es) Gea Te7s" S*transie2On ase tne 
correct interpretation, the state energy of the 
Zor l . ee 
(2s) S state can be calculated as follows: 
F 1. De Ae 
accurate energy of ls2s S;3 -2,1469 a.u. 
(28)? Is - 1s2s 1, fransition: 
279,715 em = 1.2744 a.u. 
oe. a 7 hie 
energy of (2s) S state: 0.98726 (42... 
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1 Ouse: 


Holoien (26) has set out to test this 
interpretation by constructing a Cl function for tne 
(De)7 S State. This function has an expectation 
value for the energy of -0.7985 a.u. On te basis 
of the difference between this value and the one 
calculated above, -0.8716 a.u., Holoien rejects the 
tLransition (28)7 Te Ses 2s ss aS an interpretation 
Of the corona line.” Using his Cit functions He calculates 
that other possible transition for the corona line are: 


owe 3 
gee 


(2s)° 16 - ls3p iD and (2s) ssp a2. By unl 1 cation 


Holoien claims that his CI function is a good 


Es State and that the energy 


; 2 
approximation to the (2s) 
CEE On 19S dalle ts d Close Upper bound to the trues 
energy of that state. In the following sections these 


claims are, however, rejected. 


VLG 2 Holoien's Cl function 


: Zee. ye hee es) 
In his work on the (2s) Se state (26). one’ Oo: the 


bases Holoien used to construct CI functions was basis 
No. tll listed in Table 10. As described in Chapter iV 
it was this basis which in the present work was selected 
for an initial study of the Rayleigh-Ritz variational 
method. Of particular interest here is Figure 6 in 
which the results were plotted. There is complete ag 


ment with Holoien 
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the coefficients of the CI function oe and the energy 
E, are in the present work identical to those obtained 
by Holoien. 

The main point, however, of Holoien's publication 
is the claim that with the same basis a good CI function 
is obtained for the Geers ee states, Thateftunction, 
listed as a vector of expansion coefficients, is stated 
to belong to a minimum point on an energy curve obtained 
by varying the common orbital exponent n of the basis. 
This curve is further described by Holoien as being 
identified with the (26) root of the secular equation. 
It is not clear on what basis this identification rests. 

In the present work Holoien's CI function has 
been completely reproduced but not as a function belonging 
tO a Minimum point. In fact, the function was formed 


to belong to the point marked X on the E3 curve in Figure 


' 
‘ 


Go Tha s point represents E = -0.7985 andey = 0.2375. 
Clearly that point is not a minimum of an energy curve. 
To further check for a possible minimum, calculations 
were carried out on both sides of n = 0.9375, with the 


results: 
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e detai results coincide 
Clearly the trend shown by these detailed results coincid 
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with the gross trend shown by the E, curve in the neigh- 


DOWENOOG LOf eye 0.9.37 5% 


There is no evidence of a minimum and therefore 
the present work refutes Holoien's claims for the CI 


function under discussion. As a conseguence his rejection 
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for the corona line rest on an untenable basis. 


Vis oo Ely CI method for Cie (254 I, state 


According to the variational principle the energy 
Om any thial function is an upper bound to the ground 
state energy of a system. By this same principle the 
energy of any trial function will be an upper pbownd to 
an excited state energy provided that function is oxrtho- 
gonal to the true eigen functions of all the lower 
states of the symmetry under consideration. Thus to 

x Bees? Le 

obtain the best upper bound to the energy of the (2s) Ss 
state the trial function would have to be orthogonal 
to the infinitely many lower lsns “S States. 


This orthogonality requirement can be rephrased 


in terms of vector spaces. Let the trial PuNnceLon LoT 
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in which ¢ is a symmetry adapted two-electron con- 


figuration function. Then the pr should lie 
completely outside the space of all lsns 15 state 
functions, that is outside the "lsns space", and shovld 
lie as much as possible in the space of all 2sns zs 
state functions, that is inside the "2sns space". 

in actual variational calculations this ortho- 
geonality condition can obviously not be completely 
fulfilled. However, it opens the question as to what 
eeqree partial Fulfilment cf this condition helps to 


improve an upper bound to the (2s)? ng 


State energy. 
In response to this question a CI basis was selected 
which, firstly,should span as much of the lsns space 
as possible and, secondly, should contain configuration 
zunctions lying as much as possible in the 2sns space. 

The first Fequirement is met by the 1/n basis listed 
fa table is, Since) eas shown earlier; *this is ‘the best 


basis “available at present to represent the Isns space. 


nd requirement, the success in choosing 
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the configurations of the form lsns in the j/n bas 

: ; 4 ue ' rie Gude Ont | 
obtain CI functions for 1sns 5 states suggests “ar 
analoaous choice with respect to the 2sns space, namely 
the choice of configurations of the form 2snse Thus for 
a first trial a 20 configuration basis was set up by 


5 
2 ee a. 
adding the configurations (2s), 2s3s, (3s) , (4s)", and 


104. 


ete 
o 2 
_ 


‘Be i eel’ 
be Re 


a 
ah feet. niki 


: 7 ae Ries eo Lgmot 
* id pa tei : a 
mei w 9 fue ae oan mae rf 
. Lvs ne , a Ber : Ste 4 of ‘ou we to saul ioe Fs 
tie so : tim 
Cen ene 4 OMe yy hon ee BoA gid ae ‘tebe ee aor ae 
i cM 
‘ cs 2s 7 
+) Agana eet eee, oat sicher 
a oe eee: ee at; 
R= * = ey reef, 
a ok ct - : ae A As 2) OF = #, 
! shia elbegies” SA! 2 one Tanetin ae oy 
, ie : « P ral 
eh! Sd Be. uk Pte Se | Creige Sae Sae 
a ied ad 7 4 
BE peel RAP aD PEE gpg 1 ee 2 Sea ae 
* : : f. a $ 7 : 
AMS eee es y* gayvoegs Gite 
* a J 
re 7 [e. + 
okie geile ts glee 
Dede ailatad Ooh eee 
~ ie ba | | y ; ie neg. 
} “ 
4 , ly enor au of { bs ‘5 e 
i 
‘ 7 ey biz . 
« 125504 m4) "ees ee ke, 
is t eke 
Me GS Bal ENE ak RAY peer Bd uaier seg 
he -~ 
: aT 
"y - «©. a a 
a er) : 74 , { ¢ he 4 < 4 7 © its th: 2 _ 
; army ep ’ 4 A 
¥ ae _ "1 
: ape ceog (ett Sy te shoe Bis ; oe ee 
t 


ie Te ren 43 “i 


é, 


» ve 


near f & 


a 


4 4. fi 
a ative ai 
: i] 
a.» 
ie tw 


i, 3 aera 


‘Senet? * mt ss ae 


x sncits Rass, 2¢ hie! ns a 
a * aie tp ido ae waeens 


ene aa es 


en wae 4¥iSe ie 


Gane to the 15 configurations of the i/n basis. The 
new configurations, the "2sns Pare” of the basis, were 
expressed in STO's in the same way as had been done for 
the 1/n basis, the "“lsns part" of the new basis. 

The orbital exponents of the lsns part are fixed 
as listed in Table 15. The orbital exponents of the 
2sns part of the basis were selected to serve as a 
common variational parameter to be denoted by n(2sns). 

Thus by varying n(2sns) a minimum value was sough 
for the 16th lowest eigen value, Eis: 


diagonalization of the CI matrix. This eigen value is 


obtained by the 


the first above the 15 eigen values corresponding to 

a , : yee 
the lsns “S states and is therefore expected to belong 
4. 4 1 (Bs f . : 
to the first 4sns “S state. The energy curve for Eas 
obtained by varying n(2sns) is shown in Figure 8&8. 
The minimum value for & is =-0.7338 a.u. at n(2sns) = 


0.925. The compiete set of E.'S obtained at n = 0.9 


is listed in the fourth column of Table 21 and there 


compared with the E,_'s earlier obtained from the 1/n 
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. 5 . see ES f +: 
basis (3rd column) and with the accurate energies (2nd 


column). One notes the following results: 
i) the 15 E's belonging to the lsns part of the 
ai & i 
basis (4th column) are all slightly lower than 
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as upper bounds to the accurate energies (2nd 
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t-- 
- 
— 


with E,5 = ~9.7338 the 2sns part of the basis 
does indeed establish an energy far above the 


series limit of -2.00 a.u. for the lsns 1. states. 


The question remains whether E\s can be considered 
as an upper bound to the energy of the (2s)? Is 


state, and, if so, how dependable it may be. 

Patac gsc! 1.20 Cl functions have been submitted 

to an RO analysis. The coefficients of the dominant 
configurations are listed in the second last column of 
Table 21. The last column indicates whether the dominant 
configuration is one doubly occuped orbital (by X4Xq) 

oer has two different singly occupied orbitals (by X4Xo)° 
Pea Senoced sthat stor ®9 to m4 the pattern of the..co- 


@fficients and the forms of the configurations.<is entirely 


Similar to that for the original 1/n basis as shown in 
Yable 17. Furthermore the i5 is shown having a doubly 


occupied RO for its dominant configuration. This is 
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as being an 
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to what extent the lsns part of the basis represents 

a convergence towards a complete spanning of the lsns 
Space, defined above. The better this convergence 

the more stable the lsns part of the basis can be 
considered in the sense that the removal of one lsns 
configuration should not appreciably change the lowest 
doubly excited state energy, SWORN TO test) thac 
stability, a number of calculations were performed in 
which the lsns part of the basis was reduced by omitting, 


Onewat a time;,the configurations) lel 5s) Jlsits,... is7 


uv) 


in that order. 
In Table 22 are listed the values for E({(2s°) 
btained after each reduction of the basis together with 
the coefficients of the corresponding dominant RO 
configuration. The dominant configuration retains its 


doubly occupied. character throughout the reduction of 


To show the trend of the (287) values and the 
dominant coefficients a plot is given in Figure 9. 
The curves clearly show a degree of convergence toward 
a limiting value for both E(2s") and the dominant 
coefficient. 
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Table 22 Progressive Reduction of the CI basis 


Configurations Coefficient 

Removed _ E(2s%) RO Analysis 
Ost sad Os eslehs, 
lsi5ds Oc Oa 0.9824 
1si4s —07 135.0 0.9809 
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t+islls meow OSS7A3 
isldQs -0.7419 0.9708 
-ls9s -0.7448 OQ. 9602 
‘Is8s ~0.7488 0.9603 
+1s7s Ona oo -0.9525 
+1is6s ' -=0. /639 -0.9428 
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In Table 23 some examples are given which are represen- 
tative of the overall trend throughout the reduction of 
the basis. The table has been organized to show the 
Surprising waet that, as*the basistis reduced by sone 
configuration, the eigen value removed is not the upper- 
most of the lsns part but rather the second highest. 
That uppermost 1lsns eigen value remains in the 
neighbourhood of -2 a.u., the series limit of the lsns Ig 
energies. 

Them last. column lists each Ey value reached before 


it was removed by reduction of the basis. Excepting 


the two highest E, values the largest differences 


k 
between the initial values (first column) and the final 
Veluess(last Golumn) occur an-+the third decimalss Thus 
the Ey values belonging to the lsns part of the basis 
show a remarkable stability throughout the reduction 

of the basis. 

Although further work could be done on optimization 
of the Isnsi part.ocof the basis.and on inclusion, zn the 
2Zensevacte, Of more functions to improve the radial 
andeangular correlation, the present work is adequate 
to the purpose of determining whether the (2s)- me state 
is.thesupper StaceroL, the 2797/15 om + corona line. 


The correlation energy of the ground state of He is 
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